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Abstract 

We present cross sections for production of electroweak vector boson pairs, WW, 
WZ and ZZ, in pp and pp collisions, at next-to-leading order in a s . We treat the 
leptonic decays of the bosons in the narrow-width approximation, but retain all 
spin information via decay angle correlations. We also include the effects of WW Z 
and WW j anomalous couplings. 
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1 Introduction 



At the core of the electroweak Standard Model is its invariance under the nonabelian 
gauge group SU{2) x U(l). Many aspects of this gauge structure, such as vector boson 
masses and couplings to fermions, have already been tested with high precision in a 
variety of experiments. However, the nonabelian self-interactions of vector bosons - 
in particular, triple gauge-boson couplings — are just beginning to be studied directly, 
via vector-boson pair production in e + e~ annihilation at LEP2 at CERN, and in pp 
collisions at Run I of the Fermilab Tevatron. Although thousands of W + W~ pairs have 
been collected at LEP2, they have all been produced at relatively modest values of 
the pair invariant mass, My/w ~200 GeV. On the other hand, if there are anomalous 
(non-Standard Model) vector-boson self-couplings, their effects are expected to grow 
with invariant mass, so it is useful to study vector-boson pair production at the highest 
possible energies. Vector boson pairs also provide a background for other types of physics. 
If the Higgs boson is heavy enough it will decay primarily into W + W~ and ZZ pairs 
Exotic Higgs sectors can have substantial branching ratios for charged Higgs bosons to 
decay to WZ pairs |2|. Leptonically decaying WZ pairs, W + Z — > £ + z/£' + £'~, in which 
the negatively-charged lepton is lost, form a background to a signal for strong WW 
scattering associated with the mode W + W + — > £ + vl l+ v' 0. Finally, a prime signal 
for supersymmetry at hadron colliders is the production of three charged leptons and 
missing transverse momentum [Q; a background for this process is the production of a 
W plus a (virtual) Z or 7. 

In the near future, hadron colliders will be the primary source of vector boson pairs 
with large invariant mass. Run II of the upgraded Tevatron should yield a data set 
roughly 20 times larger than that from Run I, including 100-200 leptonically decaying 
W + W~ pairs. The Large Hadron Collider (LHC) at CERN promises to increase the 
sample by another factor of 50 beyond Run II. With this increase in statistics, refined 
Standard Model predictions are essential. The leading QCD corrections (0(a s )) are 
significant (generally of order tens of percent), and hence are required to get a precise 
estimate of the overall production cross section. Also, experiments do not detect vector 
bosons, but only those leptonic decay products that fall within the experimental accep- 
tance.[] Spin correlations between vector bosons are reflected in kinematic distributions 
of leptonic momenta, which in turn influence the number of events surviving experimen- 
tal cuts. In order to properly take into account the effects of cuts on the cross section, as 
well as to study the more detailed (lepton) kinematic distributions permitted by higher 
statistics, it is important to treat the vector boson decays properly, including all spin 
correlations. 

Hadronic production of vector boson pairs in the Standard Model has already been 
studied extensively. The Born-level, or leading-order (LO) cross sections for W + W~, 
W ± Z and ZZ pair production via quark annihilation were computed twenty years ago 0. 
These cross sections were evaluated by treating the W and Z as stable particles and 
summing over their polarization states, using completeness relations to simplify the sum; 

1 Modes in which one of the vector bosons decays hadronically have been studied at the Tevatron, 
but at Standard Model levels these events are hard to separate from the QCD production of a vector 
boson plus jets 
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thus spin and decay correlations were neglected. For the spin-summed production cross 
section, the next-to-leading order (NLO, or 0(a s )) QCD corrections were obtained for 
W + W~ final states in refs. g, for W ± Z in refs. g, 0, and for ZZ in refs. [0, |]|. 

The simplest way to include the effects of vector-boson spin and decay correlations 
is to compute directly the matrix elements for the production of the four final-state fer- 
mions. In the narrow-width approximation, only 'doubly-resonant' Feynman diagrams 
have to be considered — the same class of diagrams that gives rise to the on-shell spin- 
summed cross section. Because both the outgoing fermions and the initial-state partons 
are essentially massless, and because their couplings to vector bosons are chiral, it is very 
convenient to use a helicity basis for the fermions. The tree-level helicity amplitudes 
for massive vector-boson pair production and decay into leptons were first computed in 
ref. Jl3|, which also demonstrated the significance of decay-angle correlations. This same 
approach was carried out at order a s in ref. [ H| . At 0{pc s ), there are real corrections, 



consisting of tree graphs with an additional gluon in either the initial or final state, and 
virtual corrections, consisting of one-loop amplitudes that interfere with the Born ampli- 
tude. However, the full one-loop amplitudes including leptonic decays were unavailable 
until recently, so ref. [13] included decay correlations everywhere except for the finite 
part of the virtual contribution, for which spin-summed formulae were used.0 

In this paper, we present 0(a s ) results for the hadronic production of W + W~, WZ 
and ZZ pairs, including the full lepton decay correlations in the narrow-width approx- 
imation. We rely on ref. |[l5j for all the required matrix elements, in particular the 
virtual one- loop amplitudes for qq' — > V1V2 — > 4 leptons. In order to cancel the infrared 
divergences in the phase-space integrations for the real corrections, we implement the 
general subtraction method discussed in ref. |16||. This method allows the computation 



of distributions of arbitrary (infrared-safe) observables. 

Recently, an update to vector boson pair production has been presented in ref. flT7|] . 
The corresponding Monte Carlo program, MCFM, relies on the same amplitudes |T5| and, 
therefore, also includes all spin correlations exactly to next-to- leading order in a s . MCFM 
is more complete than the program described here, in the sense that the narrow-width 
approximation is not assumed, and singly-resonant diagrams are also included. These 
additions are expected to shift the resonance-dominated di-vector boson cross sections by 
the order of several percent. Their effects are obviously much bigger in the off- resonant 
regions important for studies of Standard Model backgrounds to new physics. 

In the present paper, we first compute the di-vector boson cross sections in the 
Standard Model, both without and with a realistic set of experimental cuts. For W + W~ 
production, a jet veto is used by experimentalists to suppress backgrounds; we study the 
effect of this veto on the size of the cross section and its renormalization/factorization 
scale dependence. 

Di-boson amplitudes in the Standard Model have interesting angular dependences. 
For example, at Born level, there is an exact radiation zero in the partonic process 
qiq 2 — > W ± 'j at cos 9 = (Qi + Q2)/{Q\ — Q2), where 9 is the scattering angle of the W 

2 The virtual corrections can be divided into terms with poles in e, the parameter of dimensional 
regularization, plus residual finite terms. The pole terms have a universal form and cancel against 
infrared divergences in the real corrections, so if the real corrections include decay correlations, then 
the virtual pole terms must also, in order to get a finite answer. 
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with respect to the direction of the quark qi, and Qi^ are the quark electric charges f|. 
Similarly, it has been shown that there is a an approximate zero for qiq 2 — > W^Z at 
cos 9 = (gi + g^lidx — #2"), where g^ 2 are the left-handed couplings of the Z boson to 
the quarks fI8| . The exact tree- level zero for Wj is filled in somewhat by QCD radiative 
corrections, and also by the kinematic ambiguity associated with the undetected neutrino 
in W +, y — ► £ + i"j. Still it produces a dip in the distribution of a related variable, the 
rapidity difference y 7 — y e +, which should be visible at Run II 

Here we study the QCD corrections to the approximate WZ radiation zero. Using 
on-shell, spin-summed cross sections, ref. [10 computed the distribution in the rapidity 
difference between the W and Z bosons, Aywz = \yw — Uz\, which is a boost-invariant 
surrogate for the center-of-mass scattering angle 9. It was found that a dip in the Ay wz 
distribution persists at 0(a s ), although the dip is less pronounced than at Born level. 
Since the rapidity of the W cannot be determined on an event-by-event basis, we study 
a quantity related to Ay wz , but constructed purely out of charged lepton variables, and 
find that a dip (or at least a shoulder) is still present. However, because of the much 
lower WZ cross section, this measurement is considerably more challenging than the W^f 
case, and will probably have to wait for the LHC. 

Various types of TeV-scale physics may modify vector-boson self-interactions. With- 
out a precise knowledge of the new physics, one often parameterizes the modifications 
using anomalous coupling coefficients. We shall consider anomalous contributions to the 
W + W~Z and W + W~ r y triple gauge vertices, and their effect on various distributions 
in WW and WZ production. Similar studies have already been carried out at order a s 
and including spin correlations everywhere except in the finite virtual contributions, for 
both WW production |21J and WZ production pifl . Here we include as well the spin 



correlation effects from the finite virtual contributions. This requires matrix elements 
beyond those in ref. jTH]; however, the new matrix elements are trivial by comparison. 

The remainder of the paper is organized as follows. After outlining the computational 
techniques in section || we present results for the Standard Model production of WW, 
WZ and ZZ pairs in section |3|. We first present total cross sections for all three channels, 
without and then with a set of realistic kinematic cuts on the leptons. We consider both 
pp collisions at ^/s = 2 TeV (corresponding to Run II of the Tevatron), and pp collisions 
at y/s = 14 TeV (the LHC). We discuss the dependence of the WW cross section on a 
common renormalization and factorization scale, with and without a jet veto. For the 
WZ channel, we study the approximate radiation zero before and after QCD corrections. 

In section |] we introduce anomalous W + W~Z and W + W~j couplings, and compute 
their effect on the matrix elements in the narrow- width approximation through 0(a s ). 
We then study the effects of these couplings on a double-binned transverse energy distri- 
bution for the pair of charged leptons in W + W~ production followed by leptonic decays. 
Finally, in section [5] we present our conclusions. 



2 Computation 



It is straightforward to implement the helicity amplitudes presented in ref. and 
those including anomalous couplings (see section |), in a Monte Carlo program. The 
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tree-level and one-loop amplitudes are computed as complex numbers and the squaring, 
as well as the sum over helicity configurations, is done numerically. In order to cancel 
singularities between the real and virtual parts analytically, we use the general version 



of the subtraction method |22|] as presented in ref. [16]. Our code is flexible enough 
to compute arbitrary infrared-safe quantities, apply arbitrary cuts, and add any parton 
distribution easily^ 

In this paper we will present results for the Tevatron Run II and the LHC The former 
term refers to pp scattering at ^/s = 2 TeV, whereas the latter stands for pp scattering at 
y/s = 14 TeV. Most of our results will be presented with some 'standard cuts' which are 
defined as follows: We make a transverse momentum cut of pt > 20 GeV for all charged 
leptons. The event is required to have a minimum missing transverse momentum p™ lss , 
which is carried off by the neutrino (s). We require p™ lss > 25 GeV in the case of Vi^-pair 
production and p™ ms > 20 GeV in the case of W ± Z production. No p™ lss cut is applied 
for Z-pair production. Finally, we apply some collider-dependent rapidity cuts for the 
charged leptons. For the Tevatron we require \r]\ < 1.5, whereas for the LHC \r]\ < 2.5. 

In all results presented in this paper, we assume that the vector bosons always decay 
leptonically, i.e. the proper branching ratios of the vector boson decays into leptons 
Br(V — > //') are not included. Obviously, these branching ratios depend on which final- 
state charged leptons are included in the analysis (electrons, muons, or both). They can 
easily be added at any stage, using 

Br(Z -> e+e") = Br(Z -> fi + fT) = 3.37%, 

J2 Br(Z^z^) =20.1%, (1) 

i=e,fx,T 

Bi{W + -> e + v e ) = Br(W + -> ^v^) = 10.8%. 

These ratios implicitly incorporate QCD corrections to the hadronic decay widths of the 
W and Z. 

We use two different parton distributions, MRST(ft08a) f23fl and CTEQ(4M) [||, 
which we refer to simply as MRST and CTEQ. For both the leading and next-to-leading 
order results we shall use the same parton distributions (which have been obtained by a 
fit at next-to-leading order in a s ). The strong coupling constant is evaluated using 

. a s {Mz) ( a s (M z ) f3 1 \n(w)\ 
« S (W = 1 



W \ 71 (3q W J 

w=1 - 0o <^l iB (Ml), (2) 



7T V jJ, 

with A, = |(f C A -\N f ), fix = \CiC 2 A -{\C A +C F )N f ) ) C A = N c , C F = (iV 6 2 -l)/(2iV c ). 
The value of a s {Mz) is set equal to the value given in the respective parton distribution 
fit. Thus, we take a s (M z ) = 0.1175 for MRST and a s (M z ) = 0.116 for CTEQ. In all 
computations we have set the renormalization and factorization scales equal: fi R = fi F = 
fi. 

3 We wrote two independent programs, one in Fortran 77 and one in Fortran 90, both of which are 
available upon request. 
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The masses of the vector bosons have been set to Mz = 91.187 GeV and Mw = 
80.33 GeV. As for the coupling constants a and sin 2 #y^, we choose them in the spirit 
of the "improved Born approximation" [^5], for W pair production at LEP2. We do 
not explicitly include any QED or electroweak radiative corrections. However, we take 
into account the top-quark-enhanced corrections to the relation between Mz, Mw and 
sin 2 6 W , where the latter is defined as an effective coupling in a high-energy process, by 



using the definition [26 



. 2 _ -KajMzj 

where Gf = 1.16639 x 10~ 5 GeV -2 is the Fermi constant and the running QED 
coupling. For our numerical results we use a = a(Mz) = 1/128 and sin 2 9w = 0.230. 

The programs have been set up to allow for arbitrary values in the entries of the 
Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix that do not depend on the top 
quark. For the numerical results we take \V u d\ = \V CS \ = 0.975 and \V US \ = \V c d\ = 0.222. 

Finally we mention that we do not properly include processes where top quarks are 
involved. The amplitudes presented in ref. |H| assume massless quarks, an approxi- 
mation that is certainly not justified in the case of the top quark. The fact that we 
nevertheless include the t-channel exchange of the top quark (with \V t d\ = \V ts \ = 
and \Vtb\ = 1) for IV-pair production, therefore, results in an error. Fortunately, these 
processes are suppressed for energy scales that are not too large, either by small CKM 
matrix elements or by the small h quark distribution function. Indeed, we checked that 
the contribution of the subprocess bb — > W + W~ (treating the top as massless) is com- 
pletely negligible for Run II while it is of the order of 2% for the LHC. Furthermore, in 
the case of W ± Z production we did not include the process bg — > W~Zt. This process is 
present at next-to-leading order but is strongly suppressed by the large top quark mass, 
as well as the small b quark distribution function. 



3 Standard Model Results 



3.1 Total Cross Sections 



The total cross sections for NLO vector-boson pair production were computed long ago 
0. |, H, [T(J [U], [T2[ and have recently been updated |H] . In Tables [1] and |2] we present the 



total cross sections for the various processes at the Tevatron and LHC, for the MRST 
and CTEQ parton distributions. For the purpose of comparison we tabulated the results 
for <7 tot , the cross sections without any cuts applied, but we also give <r cut , the cross 
sections with the standard cuts defined in section At the Tevatron, the W + Z and 
W~Z total cross sections are equal by CP invariance. The cross section values are for 
the scale \i = (My 1 + My 2 )/2, where My { are the masses of the two produced vector 
bosons. Because the difference between the MRST and CTEQ results is rather small, 
and given the fact that these distributions will be updated regularly, we restrict ourselves 
in the remainder of this paper to the MRST distribution. For the same choices of input 
parameters and parton distributions, we obtain perfect agreement with the total cross- 
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ZZ 


W+W~ 


W~Z 


LO 


NLO 


LO 


NLO 


LO 


NLO 


a tot (MRST) 


1.13 


1.44 


9.52 


12.4 


1.37 


1.84 


o- tot (CTEQ) 


1.16 


1.47 


9.89 


12.8 


1.38 


1.86 


a cut (MRST) 


0.352 


0.446 


3.17 


4.22 


0.377 


0.506 


o- cut (CTEQ) 


0.362 


0.457 


3.31 


4.40 


0.385 


0.520 



Table 1: Cross sections in pb for pp collisions at ^fs = 2 TeV. The statistical errors are ±1 
within the last digit. 





ZZ 


W+W~ 


W~Z 


W + Z 


LO 


NLO 


LO 


NLO 


LO 


NLO 


LO 


NLO 


<7 tot (MRST) 


11.4 


15.2 


77.9 


115 


11.0 


19.0 


17.6 


30.1 


cr tot (CTEQ) 


11.8 


15.8 


81.3 


120 


11.4 


19.6 


18.6 


31.9 


cr cut (MRST) 


3.95 


5.31 


24.6 


40.4 


3.41 


6.44 


5.08 


9.38 


a cut (CTEQ) 


4.09 


5.51 


25.6 


42.0 


3.59 


6.72 


5.32 


9.83 



Table 2: Cross sections in pb for pp collisions at yfs = 14 TeV. The statistical errors are ±1 
within the last digit. 



sections tabulated in Tables 1 through 4 of ref. W>lu The significant differences between 
the total cross sections for the MRST distributions given in our Tables [I] and |2] and those 
in ref. |17| have their origin in different input parameters, in particular sin 2 ^^. In the 



case of the CTEQ set an additional difference is due to the fact that in ref. [T7J] the more 
recent CTEQ (5) parton distributions have been used. 

The one-loop corrections to the total cross sections are of the order of 50% of the 
leading-order term. However, as we will see below, the corrections can be much larger for 
large px or invariant mass of the vector bosons, particularly at the LHC. This is related 
to the fact that at next-to-leading order the sub-processes qg — > ViVzq have to be taken 
into account. These sub- processes generally dominate the tail of the pr-distribution (see 
e.g. ref. ||). Therefore, a scale choice like 

^ = /4 = \ipWi) + Pt(V2) + M\ + M 2 V2 ) (4) 

seems to be more appropriate. The difference between the two different scale choices 
is very small for the total cross section, since it is dominated by low pt vector bosons. 
However, for more exclusive quantities the differences can be substantial. It is therefore 
necessary to investigate the theoretical uncertainty related to the scale dependence in 
some detail. 



4 Tables 1 and 2 of ref. |17j express their good agreement with refs. [|[ |T(j, [T2| for the older HMRSB 
parton distributions; Tables 3 and 4 are for the Tevatron Run II and LHC with the MRST and CTEQ (5) 
sets. 
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Figure 1: Scale dependence of <r cut , the cross section for VT-pair production at 
the Tevatron with standard cuts. The scale is given in units of Mjy. We show 
the LO, NLO and NLO with jet veto curves. The inset shows the three curves 
normalized to 1 at /i = M\y. 



To start with, in Fig. |l] we consider the scale dependence of the cross section for Im- 
pair production at the Tevatron. We apply our standard cuts and vary the scale around 
/i = Mw The leading-order scale dependence is entirely due to the decrease of quark 
distribution functions q(x) with increasing factorization scale for moderate x. The NLO 
result does indeed have a reduced scale dependence. We also show the NLO result with 
an additional cut on the transverse hadronic energy, _E^ ad < 40 GeV, implemented at the 
parton level. In the remainder of this paper we refer to this additional cut as jet veto, 
even though it does not exactly correspond to a jet veto applied by experimentalists. This 
additional cut reduces the scale dependence further. These results, together with the 
modest one-loop corrections we find below for kinematic distributions at the Tevatron, 
lead to a rather satisfactory description of VT-pair production at this collider. 

The situation is somewhat more delicate for the LHC. The one-loop corrections can 
be huge in the tails of the distributions. We therefore investigate the scale dependence 
in a more detailed way. We again consider the cross section with and without the same 
jet veto, E^? d < 40 GeV. We also consider the cross section with an additional pair of 
cuts on the transverse momenta of the charged leptons. We require the larger of the 
two transverse momenta of the charged leptons, p™ ax , to be bigger than 200 GeV and 
the smaller, p™ in , to be bigger than 100 GeV. As discussed above, it is therefore more 
appropriate to vary the scale around /i as given in eq. (^) instead of /i = Mw- 

The purpose of these additional cuts is to investigate the scale dependence of the 
cross section in the region where larger corrections are expected. Indeed, the one-loop 
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Figure 2: Scale dependence of <7 cut , the cross section for VF-pair production 
at the LHC with standard cuts. The scale is given in units of // s t as defined in 
eq. (|j). We show the LO, NLO and NLO with jet veto curves without additional 
cuts (left) and with an additional cut p™ x (£) > 200 GeV and p™ in (£) > 100 GeV 
(right). The insets show the curves normalized to 1 at p, = p st . 



correction to the LHC cross section for fi = /i st increases from 60% to 80% as the 
additional set of cuts is applied. Before applying the additional cuts, the situation is 
very similar to the Tevatron. The scale dependence at LO is reduced at NLO, and is 
reduced even further if the jet veto is applied. For the high px case this is not quite 
true. The leading-order result is surprisingly stable under scale variations. This feature is 
somewhat artificial, however; we have checked that it does not hold if the additional cuts 
are changed to e.g. p™ ax (f) > 400 GeV and p™ in (£) > 200 GeV. In this case, the leading- 
order cross section decreases with increasing scale. With the cuts p™ ax (£) > 200 GeV 
and p T an {£) > 100 GeV we just happen to be close to the transition from a rising to 
a falling leading-order cross section, a transition which is associated with the different 
behavior under evolution of the quark distribution functions at small x vs. moderate x. 

The reduction of the scale dependence of the next-to-leading order results when a 
jet veto is applied seems to be quite general though. This situation is a bit paradoxical 
because the cross section with a jet veto is less inclusive and, therefore, expected to 
be more sensitive to large logarithms created by incomplete cancellation of the infrared 
singularities. On the other hand, any subprocess appearing at NLO in pp —>■ WW 
produces additional hadronic energy in the final state. Thus, a cut on -E^ad naturally 
suppresses the one-loop corrections and tends to stabilize the perturbative expansion. 
This effect competes against the stronger sensitivity to large logarithms. Apparently, for 
the jet veto we applied, E^ ad < 40 GeV, the subpro cess-suppression effect still dominates. 

3.2 W + W~ 

In this subsection we present some results for kinematic distributions for the processes 
pp —* W~W + — > £~u£' + u' and pp — > W~W + — > £^v£' + v' . Similar studies have been 
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carried out earlier (see e.g. ref. fL4|)- Throughout, we apply our standard cuts as defined 
in section |2|. As an illustration, we have chosen eight variables, four p^-like quantities 
and four angular distributions. They all are defined in terms of observable momenta. 
The pr-hke distributions are defined as follows: 

transverse momentum of the negatively charged lepton. 

Ma'- invariant mass of the lepton pair. 



pmiss. m i ss i n g transverse momentum, J (pr(£~) + Pt(£' + ) +Pr(jet)) 2 . 

^max. m aximal transverse momentum of the two charged leptons, max{pjn(£~) , pj-(£' + )} . 

In the case of the Tevatron, the Pt{^~) and Pt{^' + ) distributions are equal. However, 
this is not true for the LHC. 

The four angular distributions we considered are defined as follows: 

rj(£~): rapidity of the negatively charged lepton. 

Ar](£): rapidity difference between the leptons, 1](£~) —i](£' + ). 

cos6*«: angle between the leptons, cos(l[p(£~) , p(£' + )}) . 

cos0«: transverse angle between the leptons, cos(l\pr(£~), pt{£' + )])- 

For massless leptons, the true rapidity y(£) = | ln ^^|^|^ is equal to the pseudora- 
pidity r](£) = — ln(tan(#/2)), so we refer to both as rapidity. For the Tevatron, we have 
to specify the proton direction: it corresponds to 9 = 0. 

At the Tevatron, the rapidity distribution for £' + can be obtained trivially from that 
for £~ by changing the sign of rj. For the LHC, there is no such simple relation between 
the two distributions. Note also that the Ai](£) distribution is not symmetric around 
rj = for the Tevatron; it is symmetric for the LHC. The cos 8 a observable has been 
investigated in ref. p7[ in the context of Higgs boson detection in the intermediate Higgs 



mass range tuh = 155-180 GeV. It is therefore particularly interesting to consider the 
effect of the QCD corrections to this observable. 

In Figs. ^| and £| we show the distributions for the Tevatron, evaluated at n = M w 
and with our standard cuts applied. Recall that the branching ratios for the leptonic 
decay of the vector bosons are not included. The perturbative expansion for angular-type 
distributions is generally better behaved than that for steeply falling p^-like distributions. 
The latter often suffer from large NLO corrections in the tails of the distribution. The 
insets of Fig. |3] show the ratio da NLO /da LO . They indicate that for the Tevatron, the 
NLO corrections are not too large, except for the p™ lss distribution. 

The large corrections to the p™ lss distribution reflect a suppression of the Born-level 



cross section when p™ lss is large |2(| : The only way to get a large p™ lss at Born level is 
to have a large WW invariant mass. Then the two neutrinos are almost back-to-back 
in the WW rest frame. Also, for large invariant mass the Standard Model WW helicity 
amplitudes are dominantly those where the W + and W~ have opposite helicities. The 
V — A decay of the W bosons then implies that two neutrinos tend to carry roughly the 
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Figure 3: Differential cross sections in pb/GeV for W~W + production at the 
Tevatron for the pr-hke variables pr(£~), M«, p™ lss and p™ ax defined in the text 
at LO (dashed curves) and NLO (solid curves), with fi = My/. Standard cuts 
have been applied and the branching ratios for the leptonic decay of the vector 
bosons are not included. The insets show the ratio da NLO /da L °. The units on 
the horizontal axes are GeV. 
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Figure 4: Differential cross sections in pb for the Tevatron for the W~W + angu- 
lar variables defined in the text at LO (dashed curves) and NLO (solid curves), 
with \i = Mm/. Standard cuts have been applied and the leptonic branching 
ratios are not included. 
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same fraction of the W momentum, i.e. they have cancelling transverse momenta. NLO 
QCD corrections have a big effect at large p™ lss because they allow a recoiling final-state 
parton to spoil the px balance of the two neutrinos. The presence of anomalous couplings 
can also have a big effect in this region, by relaxing the helicity anti-correlation of the 
two W bosons @. 

In Figs. H and [| the same two sets of distributions are displayed for the LHC, again 
for = Mw The NLO corrections to the angular distributions are again modest. The 
Pr-like distributions, however, can have much larger corrections, particularly in their 
tails, where the NLO result can easily exceed the leading order result by a factor of 5. 
The p™ lss distribution has the largest corrections of all, for the reason mentioned above, 
and the factor can be 20 or more. As mentioned before, the huge corrections in the tail 
of the generic p^-like distributions have to do with the subprocess qg — > WWq which 
dominates in this kinematical region. It may be argued that a scale choice as in eq. (B) 
is mandatory in this case. However, we have checked that such a scale choice does not 
lead to a substantial improvement. 

The problem is that in a NLO calculation for pp — ► WW the partonic process qq — ► 
WW is included at NLO, but the qg — > WWq subprocess is included only at LO. 
Therefore, in a kinematical region where the latter dominates, the calculation presented 
in this paper is effectively only a leading order calculation. The only fully satisfactory way 
to improve the theoretical prediction in such cases is to include the one-loop corrections 
to the subprocess with a qg in the initial state. These 0(a 2 s ) corrections correspond 
to a NNLO contribution to pp — > WW. At the same order in a s , there are two new 
subprocesses: gg — > WW at one loop, and gg — > WWqq at tree level. Due to the large 
gluon density at small x, these contributions are also expected to be important for LHC 
energies, and would have to be included if a reliable prediction for the tails of the p^-like 
distributions were required. 

As discussed above, a way to suppress the partonic processes with gluons in the 
initial state is to impose a jet veto. The effect of the cut E^ d > 40 GeV on the py-like 
distributions is even more dramatic than its effect on the total cross section. This can 
be seen in Fig. |5], where we also show the NLO curves with the jet veto (dot-dashed 
lines). Indeed, the one-loop corrections are very small throughout. We conclude that if 
a reliable theoretical prediction is desired for the tail of a p^-like distribution, and only 
a next-to-leading order program is available, then a jet veto is unavoidable. 



3.3 W ± Z 



In this subsection we study WZ production, followed by leptonic decay of each boson. In 
particular, we examine the effects of QCD corrections on the approximate radiation zero 
at cos6> = (gi+g^)/ (91—92) HH- Since the precise flight direction of the W boson is not 
known, due to the uncertainty in the longitudinal momentum carried by the neutrino, we 
simply choose to plot a distribution in the (true) rapidity difference between the Z boson 
and the charged lepton coming from the decay of the W, Ayze = yz — He- This quantity 
is similar to the rapidity difference Aywz = \yw ~ Uz\ studied in ref. 
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but uses only 

the observable charged-lepton variables. It is the direct analog of the variable y 7 — ye+ 



considered in ref. [[Tj| for the case of W'y production. (It is possible to determine cos 9 
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Figure 5: Differential cross sections in pb/GeV for W~W + production at the 
LHC, for the py-like variables defined in the text at LO (dashed curves) and NLO 
(solid curves), with fj, = My/. Standard cuts have been applied and the leptonic 
branching ratios are not included. Also shown (as dot-dashed lines) are the NLO 
curves with a jet veto Ej^ < 40 GeV. The insets show the ratio da NLO /da LO . 
The units on the horizontal axes are GeV. 
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Figure 6: Differential cross sections in pb for the LHC for the W~W + angular 
variables defined in the text at LO (dashed curves) and NLO (solid curves), with 
fj, = M\y. Standard cuts have been applied and the leptonic branching ratios are 
not included. 
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Figure 7: For WZ production at Run II, followed by leptonic decays of both 
the W and Z bosons, we plot the distribution, in picobarns, in the rapidity 
difference between the Z and the charged lepton t from the decay of the W, 
^Vzi = Uz — Ve- Leptonic branching ratios are not included and the scale has 
been set to \i = (Mw + Mz)/2. The basic cuts used are pr(€) > 20 GeV and 
| r){t)\ < 2 for all three charged leptons, and a missing transverse momentum cut 
of p™ lss > 20 GeV. We plot the Ayze distribution with these cuts, and also after 
imposing an additional cut on the W decay lepton, E{£) > 100 GeV; the latter 
curves have been scaled up by a factor of 5. The dashed curves are Born-level 
results; the solid curves include the 0(a s ) corrections. 



in the W'j or WZ rest frame, by solving for the neutrino longitudinal momentum using 
the W mass as a constraint, up to a two-fold discrete ambiguity for each event [EH]. 



However, ref. |19| found that the ambiguity degrades the radiation zero — at least if 



each solution is given a weight of 50% — so that the rapidity difference y 7 — yi+ is more 
discriminating than cos#.) 

As one can see from Fig. ^, there is a residual dip in the Ay Z e distribution, even at 
order a s . This dip can easily be enhanced by requiring a minimal energy for the decay 
lepton from the W . In Fig. [5], we have chosen E(£) > 100 GeV. Note that these curves 
are scaled up by a factor 5. Unfortunately, only a few tens of WZ — > leptons events 
are expected at Run II of the Tevatron, so the observation of such a dip will be rather 
difficult prior to the LHC. 



4 Anomalous VF+VF Z and W + W 7 Couplings 
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4.1 Triple Gauge Boson Vertices 

New physics may modify the self-interactions of vector bosons, in particular the triple 
gauge boson vertices. If the new physics occurs at an energy scale well above that being 
probed experimentally, it can be integrated out, and the result expressed as a set of 
anomalous (non-Standard Model) interaction vertices. 

Here we consider anomalous W + W~Z and W + W _ 7 trilinear couplings, and their 
effects on the hadronic production of WW and WZ pairs up through order a s . The 
most general set of Lagrangian terms for WWV, V G {Z, 7}, that conserves C and P 



separately, is (see e.g. J29|, |30| ) 

C/gwwv = ig\{W; v W^ - W, U W* »V V ) + m v w;W v V^ + i^-W^W^ , (5) 

where X^ u = d^X v — d v X^ and the overall coupling constants gwwv are given by gww-y — 
— e and gwwz = — e cot 6 W respectively, with 9 W the weak mixing angle. The Standard 
Model triple gauge boson vertices are recovered by letting g\ — > 1, k v — > 1 and \ v — > 0. 
The coupling factors can be written in terms of their deviation from Standard Model 
values: g\ = 1 + AgY and k v = 1 + Ak . Electromagnetic gauge invariance requires 
gl = 1, or AgJ = 0. Sometimes other constraints are imposed on the couplings. For 
example, if one requires the existence of an effective Lagrangian with SU(2) x U(l) 
invariance, and neglects operators with dimension eight or higher, then the number of 
independent coefficients in eq. © is reduced from five to three, 

A Z a W<j> X7 X Z 

A 9i = — TE~ ' A 7 = A = a W} 

COS z By/ 

a 7 a z sin 2 6^ , . 

cos z 0w 

where aw, &W4> an d o.B4> are coefficients of the dimension-six operators in this effective 



Lagrangian p0| . If one arbitrarily supposes that aw 4, = cts^, one arrives at the so-called 



HISZ scenario [31], with only two independent couplings. 

The momentum-space vertex W~ '(q)Wp \q)V^{p) (where p + q + q = 0) corresponding 
to eq. (|^) can be written as 

T^(q,q,p)/g WWV = <fg^(g\ + n v + \ V -£r) - q^(g? + K V + \ V 



+ (q" ~ q^ 



M W J - 1 * V 1 M w 



(7) 



Here we have used momentum conservation, and the fact that the terms q a ,q@ and p^ 
can be neglected, but we have not imposed on-shell conditions on the vector bosons. As 
it stands, this vertex will eventually lead to a violation of unitarity. To avoid this, the 
deviations from the Standard Model, Ag\ , Ak v and A , have to be supplemented with 
form factors. Since the form factors are supposed to be produced by unknown physics, 
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the form they should take is a priori somewhat arbitrary. We choose a conventional 
dipole form factor, i.e. 

Aj v , M Ar v , A ^ V A v ; XV (8 ) 

91 (1 + s/A 2 ) 2 ' (1 + s/A 2 ) 2 ' (1 + s/A 2 ) 2 ' lj 

where s is the invariant mass of the vector boson pair, and A is in the TeV range. 



4.2 Tree, Virtual and Bremsstrahlung Amplitudes for WW 



Replacing the Standard Model vector-boson-vertex by the more general vertex given 
above results in the some modifications of the primitive amplitudes presented in ref. . 
We use the same notation in this paper and refer the reader to ref. for more details. 
The box-parent primitive amplitudes are not affected by changes in the trilinear vector- 
boson-vertex. The change in the triangle-parent primitive amplitude can be obtained 
by simply computing the one tree-level diagram with the new vertex. Since the vertex 
is no longer symmetric in the exchange W V, we get slightly different results for the 
different final states. For the WW final state, the new tree amplitude, which replaces 
A tree,b in eq _ ^.9) of ref. fL|, is 



.4 



trcc/3 



2 Sl2 S34 S56 



{g\ + k v + A y )«13) [24] (6|(1 + 2)|5) + (16) [25] (3|(5 + 6)|4» 



A 



v 



(1|(3 + 4)|2)( 2g\ (36) [45] + _ (3|(1 + 2)|5)(6|(1 + 2)|4) 



(9) 



For the limit g\ -»• 1, n v -»• 1, A T 



0, we recover the Standard Model result, 



.4 



tree, b 



Sl2 s 34 ^56 



- (36) [45] (1|(5 + 6)|2) + (13) [24] (6](1 + 2)|5) 
+ (16) [25] (3|(5 + 6)|4)" 



Sl2 ^34 S56 



(13) [25] (6|(2 + 5)|4) + [24](16) (3|(1 + 6)|5) 



(10) 



With the above results we also get immediately the one-loop primitive amplitude 
for anomalous couplings. Assuming the usual decomposition into finite and divergent 
pieces, the finite pieces are vanishing, as in the Standard Model case, and the divergent 
pieces are still given by crA tree ' B V where 



V 





7 
2 



(11) 



The result for the bremsstrahlung diagrams with an additional positive helicity gluon 
radiated off the quark line is given by 



tree,B 



2 (17) (72)s 3 4S 56 ti27 
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{ffX + * v + A l/ )((13) (1|(2 + 7)|4)(6|(3 + 4)|5) - (16) (1|(2 + 7)|5)(3|(5 + 6)|4» 



+ (1 1 (3 + 4) (2 + 7) 1 1) ( 2g\ (36) [45] + (3| (4 + 6) 1 5) (6| (3 + 5) |4) 



M 2 W 



(12) 



replacing A t j ee ' b in eq. (2.22) of ref. [I5|. The result for a negative helicity gluon can be 
obtained by the usual flip operation fll5" |. 

We also have to modify the prescription in ref. |T5j] for dressing the above primitive 



amplitudes with electroweak couplings. Since g\ , k v and X v are relative couplings, i.e. 
the overall coupling gwwv has not been changed, the dressing with electroweak factors 
is almost identical to the Standard Model case. The only subtlety is that both Z and 7 
appear as intermediate states in WW production. In the coefficient functions 

s i2 (l =F 2Qsin 2 9 W ) 



C Lr} = ±2Qsin z 9 w + 



S12 - Mf 



(13) 



C Rr} = ±2Q sin 2 6 W =f 2Q sin 2 



8\2 



Sl2 



M 2 



defined in ref. [15], the first term (±2Q sin 2 8w) is from the intermediate 7, while the 
second term is from the intermediate Z. Correspondingly, we should set V = 7 (V = Z) 
in eqs. @ and (|i~2l), when they are dressed with the first (second) term in C,i,,«j. 
Otherwise, all the prefactors remain the same and only the 'new' primitive amplitudes 
have to be plugged in. 



4.3 Tree, Virtual and Bremsstrahlung Amplitudes for WZ 

For the WZ final state, the new tree primitive amplitude is 



.4 



trcc,B 



2 Sl2 S34 S56 



(9? + k z + A z 512 



M 2 W 



(36) [45] (1|(5 + 6)|2) 



+ {g z + K z + X z )(W) [25] (3|(1 + 2)|4) 



+(6|(3 + 4)|5) (2gf (13) [24] 



Ml 



(3|(5 + 6)|2)(1|(5 + 6)|4) 



The new bremsstrahlung primitive amplitude is 



—I 



2 (17) (72) 534*56*127 



gf + n z + y 



M 2 W 



(36) [45] (l|(5 + 6)(2 + 7)|l) 



+ {gf + k z + \ z ) (16) (1|(2 + 7)|5)(3|(5 + 6)|4) 
+ (6|(3 + 4)|5)(-2^f (13) (1|(2 + 7)|4) 
\ z 



M 2 W 



(3|(5 + 6)(2 + 7)|l)(l|(5 + 6)|4) 



(14) 



(15) 



18 



TTTnax \ pimin 
£L/j< \ JZ/^p 


20 -38.1 


38.1 - 72.5 


72.5 - 138 


138 - 263 


263 - 500 


20 - 38.1 
[Pb] 


1.07 

1.07 1.10 
1.03 










38.1 - 72.5 
[Pb] 


1.62 

1.62 1.61 
1.54 


0.77 

0.76 0.77 
0.74 








72.5 - 138 
[lO^pb] 


1.67 

1.60 1.76 
1.49 


3.30 

3.31 3.56 
3.15 


1.30 

1.32 1.50 
1.27 






138 - 263 
[10 2 pb\ 


0.63 

0.57 2.43 
0.49 


1.27 

1.20 4.18 
1.09 


3.07 

3.16 6.90 
2.98 


1.08 

1.14 2.32 
1.09 




263 - 500 
[10- 4 pb] 


0.7 

0.6 22 
0.5 


1.2 

1.1 37 
0.9 


2.6 

2.5 55 
2.2 


9.0 

9.8 63 
9.0 


2.3 

2.5 9.5 
2.4 



Table 3: Double-binned Et cross sections for pp — > W + W~ — > leptons at ^/s = 2 TeV. The 
three numbers in the left column in each entry are the Standard Model results for the scales 
/i = fi s t x {1/2,1,2}, with jUst given in eq. (|j). The fourth number in the entry is the cross 
section with anomalous couplings, as defined in the text. 

In this case, since there is only an intermediate W, the dressing with electroweak coupling 
factors is indeed identical to the Standard Model case. 



4.4 Numerical Results 

More systematic studies of the effects of anomalous couplings on hadronic production of 



vector boson pairs have been carried out elsewhere |32|, [33|, [34], Here we merely 

consider one sample distribution. 

The effect of anomalous couplings is enhanced for gauge bosons that are produced 
at large transverse momentum. In order to exploit this feature, the DO collaboration 
considered a double-binned Et spectrum for the charged leptons coming from the decay 
of W pairs ]35| . Any deviation from the Standard Model should be more pronounced in 
the high Et bins. 

We have computed a similar double-binned Et spectrum at NLO for Run II of the 
Tevatron. As in ref. |35j we compute for each event the larger and smaller transverse 
energies of the two leptons, Z?™ ax and _E™ m . (E T is equivalent to Pt(^), of course.) We 
impose our standard event cuts and then bin each Et into five bins with the following 
limits in GeV: 

E T = {20, 38.1, 72.5, 138, 263, 500}. (16) 

In order to get a feeling for the theoretical uncertainties, we repeated the Standard 
Model computation for three scales, /i st /2, /i st and 2/i st , where /i st is given in eq. (Hf). We 
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Figure 8: Cross section a [fb/10] for P^-pair production at the Tevatron with 
no anomalous couplings (left) and for the HISZ scenario with aw = c*W4> = 
OiBcf> = 0.1 (right). Standard cuts have been applied and the scale has been set 
to \i = fj, s t as defined in eq. (||). 



computed the same NLO double-binned cross section including the effects of anomalous 
couplings. As an illustration we have chosen the HISZ scenario with aw = ctw4> = <^B<j> = 
0.1. This corresponds to the following values for the anomalous couplings appearing in 
the Lagrangian in eq. (|5|): 

Ag"{ = 0; Ag z = 0.13; X 7 = \ z = 0.1; Ak 7 = 0.2; Ak z = 0.07. (17) 

The form factors have been chosen according to eq. ([8]) with A = 2 TeV. At present 
these values are still consistent with LEP2 bounds |36|]. Table [3] presents the results. 
As usual, the leptonic branching ratios are not included, and we use the MRST parton 
distributions. The three numbers in the left column of each box give the Standard Model 
result for the three choices of /i, with [i increasing from top to bottom. The fourth number 
in each box is the result with the anomalous couplings included and /i = /i st . Units for 
the cross sections for each row are given in square brackets. As expected, the results for 
the Standard Model and anomalous cross section are very similar for the low Et bins. 
For the high Et bins, however, the differences are large and certainly much bigger than 
the most conservative estimate of the theoretical error. 

The same results are also shown in Fig. || where we plot the natural logarithm of the 
total cross section (in units of 10 _1 fb) for each bin for the Standard Model with /x = /i st 
(left) and the HISZ scenario with \i = // st (right). Again, the significant differences in 
the high Et bins become apparent. 
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5 Conclusions 



We have presented a general purpose Monte Carlo program that is able to compute 
any infrared- safe quantity in vector boson pair production at hadron colliders at next- 
to- leading order in the strong coupling constant a s . This program generalizes previous 
calculations (with the exception of the recent ref. in that the spin correlations are 
fully taken into account. The decay of the vector bosons into leptons was included 



in the narrow-width approximation, whereas MCFM [pj]] also includes the singly-resonant 
diagrams needed to go beyond this approximation. For the total cross-section, computed 
in the narrow-width approximation, our results agree perfectly with those of MCFM for 
the same choice of input parameters. 

As an illustration of the usefulness of the program we presented several distributions 
for the Run II at the Tevatron and the LHC However, we refrained from performing a 
detailed phenomenological analysis; this is probably best done once the data are available. 

In addition to Standard Model processes, we considered also the inclusion of anoma- 
lous couplings between the vector bosons. We presented the one-loop amplitudes for a 
generalized trilinear W + W~Z and W + W~ , y vertex. The inclusion of these amplitudes 
into our program allows the calculation of anomalous effects at next-to-leading order, 
with full spin correlations. These effects are shown to be very prominent for large trans- 



verse momentum of the gauge bosons, confirming results of ref. p0| . Such an analysis at 



Run II of the Tevatron should yield improved bounds on anomalous couplings, although 



for large improvements one probably has to wait for the LHC [pO 
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